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ABSTRACT 

Simultaneous measurements of sea surface elevation and horizontal 
and vertical particle velocities at 39 and 69 cm elevations in the column 
of water of 130 cm total depth were made inside the surf zone. Also, 
the offshore sea surface elevation at this location was measured for 
purposes of comparison. The velocities were measured using electro- 
magnetic flow meters, and the sea surface elevation was measured using 

pressure wave gauges. Probability density functions, pdf, were determined 

$ 

for each record. The pdf's for the sea surface elevation and particle 
velocities inside the surf zone were highly skewed. Spectral computations 
show that the range of significant energy was between 0.05 and 0.6 hertz. 
The phase angle was compared to linear wave theory and shows a shifting 
of phase for the horizontal velocity with sea surface elevation from 0 
degree at low frequency to 90 degrees at higher frequencies. The energy- 
density spectra show that the horizontal component is approximately 95% 
of the total kinetic energy of the surf zone. In the range of significant 
energy, a coherence of about 0.9 was found for the sea surface elevation 
and particle velocities which indicates that the particle motion inside 
the surf zone is for the most part wave-induced. 
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I. INTRODUCTION 



A. REVIEW OF PREVIOUS WORKS 

With reference to the surf zone only the more general features 
and characteristics of the motion of the fluid are at present understood. 
This is due to the lack of understanding of waves after the breaking 
point. Most of the knowledge at this time is based on approximations 
of wave theory or empirical relationships . 

The difficulties are practical as well as theoretical. The surf 
zone is a very hostile environment to work in and a very difficult one 
to reproduce properly in the laboratory. Another major difficulty has 
been a lack of instrumentation to make velocity measurements for wave- 
induced motion and turbulence. Only a few experiments have been con- 
ducted to measure and study the kinematics of water particle motions 
due to waves in shallow and intermediate water, and almost none for 
the kinematics of water particle motion inside the surf zone itself. 

Among the laboratory experiments the most notable is a study by 
Iversen (1953) in which photographic techniques were used to obtain a 
Lagrangian description of water particle motion. However, due to the 
slope of his model beach, spilling breakers were not considered. 

Field measurements have been made by a number of authors using 
a variety of instruments. Inman (1956) measured the drag force to infer 
water particle motion. Walker (1969) used propeller type flow meters. 
Miller and Zeigler (1964) used meters based on acoustic principles; 
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they compared their measurements with higher order wave theory and 
found some qualitative agreement. Thornton (1969) used an electro- 
magnetic flow meter and presented his results in the form of spectra. 

His results showed a gradual decay of energy across the surf zone and 
a shifting of energy to the higher frequency turbulent region of the 
spectrum. 

Some attempts have been made to theoretically describe the 
kinematics of the surf zone. One of them by Collins (1970) describes 
the probability distribution functions of the wave characteristics of 
wave height and period for the region inside the surf utilizing the hydro- 
dynamic relationships for shoaling and refraction. 

B . OBJECTIVE 

The objective of this research was to make preliminary studies on 
the kinematics of the water particle motion within the surf zone and 
within breaking waves. With this purpose in mind, simultaneous measure- 
ments were made of the instantaneous sea surface elevation and of 
horizontal and vertical particle velocities at different elevations in the 
same column of water in the surf zone, and of the offshore sea surface 
elevation. The probability density functions and spectra of the wave 
and particle velocity measurements were determined. 
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II. NATURE OF THE PROBLEM 



A. CHARACTERISTICS OF THE SURF ZONE 

Theories developed to describe characteristics of waves, such as 
wave height, period, and particle velocity, can be applied reasonably 
well to deep water waves, and with less accuracy to shoaling waves up 
to the point of near-breaking conditions. However, upon breaking the 
waves lose their ordered character and can no longer be described 
analytically. Each type of breaker, spilling, plunging, or surging, 
(depending upon the beach slope and deep water wave steepness) causes 
a different type of flow and consequently a different distribution of 
energy in the surf zone. Furthermore, moving boundaries at the bottom 
and surface of the sea together with the unsteady flow motion, make the 
surf zone a very difficult place for the application of a theoretical model. 

However, results of experiments by Thornton (1969) suggest that 
the kinematics is not as disorganized as one might be led to believe. 

In order to develop a theoretical model it is necessary to have at least 
some a priori knowledge of the characteristics of the wave induced 
motion and the structure of the turbulence. 

If it is desired to study the details of surf zone flow, it is 
necessary to have a very large model or to take measurements in situ 
with appropriate instruments for measuring instantaneous velocities. 

A limitation inherent in laboratory studies is space. The most commonly 
occurring breaking wave is of the spilling type which is observed on 
very flat beaches; this calls for very long wave tanks. 



9 



Direct field measurements can overcome these problems. A 
difficulty is that the large scale conditions peculiar only to the place 
where the data is taken may affect the mean flow in such a way that it 
is applicable only to a particular set of conditions. However, small 
scale features, such as turbulent motion in the breaking waves, are not 
affected directly by the large scale conditions and their study produces 
better results when done in the field. Furthermore, the process generally 
can be assumed to be stationary, that is, the mean and variance do not 
change with time. This is a good assumption for short measurements on 
the order of 2 0 minutes as done in this research. 

In general, the surface waves and wave induced water particle 
motion have the characteristics of random phenomena, therefore it is 
necessary to utilize statistical procedures to describe them. The basic 
analyses and their application to the present problem are described in 
the following sections . 

B . STATISTICAL ANALYSIS 

1 . Probability Density Function 

The probability density function for random data describes 
the probability that a particular process will assume a value within some 
defined range at any instant of time. The probability that the sample 
time history record x(t) assumes a value between x and (x + Ax) will 
approach an exact probability description as the length of the sample 
record, T, approaches infinity, and as Ax approaches zero. In equation 
form: 
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p(x) 



AX »- 



Lim 



Lim 



1 

T 




where is the total amount of time that x(t) falls inside the range 
(x, X + Ax), and T is the total observation time. The probability 
density function is always a real, non-negative function. 

2 . Energy-Density Spectrum 



transform of the auto -covariance function of the record time series. The 
auto -covariance function describes the general dependence of the values 
of the data at one time on the values at another time. It is expressed as 
a function of lag time r as 



The quantity ( r ) is always a real-valued even -function, 
with a maximum at r = 0. 

A Parzen lag window is applied to the covariance functions to 
account for the finite length of the record in computing the spectra. The 
Parzen lag window has the advantage of no negative side lobes and 
maintains the sample coherence between its theoretical values of 0 and 
1 . This eliminates any numerical instability problem that may arise in 
the analysis procedure, particularly when calculating cross-spectra . 

The Parzen lag window is given in Bendat and Piersol (1966) by the 
formula 



The energy-density spectrum is computed from the Fourier 




t+T 



Xj^(t) x^ (t+ r ) dt 



11 



P(r) 



^9 r ? 

1 - 6 {—) + 6 {—r 

m m 



2 

0 



-I 3. 



1 “ ( X ) 
m 





T <■ T 

m 



T 



> T 



m 



where T is the total time lag. 
m ^ 

The Fourier transform of the modified auto-covariance function is 
the energy -density spectrum, given by 



00 

CDii(f)= j Pi ^ ^ dr 

-00 



This spectrum describes the partitioning of energy-density by 
frequencies. The area under the wave height spectrum is proportional 
to the total potential energy -density of the sea surface elevation. The 
area under the curve of the velocity spectrum is proportional to the 
kinematic energy -density . In both cases the total area under the 
spectrum is equal to the mean square value of the random variable and 
therefore equal to the variance. 

3 . Cross-Spectral Density Function 

The cross-spectral density is defined as the Fourier transform 
of the cross -covariance function. The cross -covariance function for 
two sets of random data describes the general dependence of the values 
of one set of data on the other. It is similar to the auto-covariance 
function except that two records are lagged against each other instead 
of a single record lagged against itself. The cross-covariance function 
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tends to the exact value as T approaches infinity. In equation form 

t + T 

Xj^(t) (t + T ) dt 

The function ^ ^ always a real-valued function, but does not 

necessarily have a maximum of r = 0, or is an even-function, as was 
the auto -covariance function. 

The cross -spectral density function is expressed as 

r * 

Oi 2 (f)= / P( )<^12^ r )e"'^ dr 

-00 

where P( r ) is the Parzen function defined above. 

Because the cross -correlation function is not an even-function, 
the cross -spectrum is generally complex and can be defined in terms of 
its real and imaginary parts 




(Diz^f) 
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where the real part is called the co-spectrum 
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P( r ) ^ "’’^12^" (2 3T f T ) dr 



and the imaginary part is called the quadrature -spectrum 
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4. 



Phase Angle 



The cross-spectral density function can be expressed in 
complex polar notation as 



Ol2«) = 






-J 



^ 12 ® 



where the phase angle, C ^ 2 ^^^ ' average angular difference by 

which the cross -correlated components of x (t) lead those of x (t) in 

2 1 

each spectral frequency band. The phase angle is calculated using the 
expression 



C 



-1 

<= 12 ® 



where C „ (f) and Q,«(f) are the cross -spectral components described 

±4^ X 4^ 



above. The phase angle is bounded by the values 



— ir < C ]^2 ® ^ 



5. 



Coherence Function 



The coherence function gives an expression for the correlation 
of two random variables as a function of frequency. It is defined by 



the equation 



Ti2 (0 = 



<Di 2 (f) 



Oil® O22® 



and has a range of values 

0 < T^2® ^ ^ 

2 

In the ideal case of two completely coherent records, ^12® 

maximum value of unity. As the correlation between the records decreases 
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the value of T decreases, and reaches the minimum value of 
zero when the records are completely incoherent or statistically 
independent. 
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III. INSTRUMENTATION 



A. FLOW METER 

The water particle velocities were measured using two Engineering- 
Physics Company water current meters, type 6130, see Figure 1. The 
current meter operation is based on the electromagnetic induction principle. 
A symmetrical magnetic field is generated in the water by a driving coil 
imbedded in the probe. When water in the vicinity of the probe has a 
velocity relative to the magnetic field, an electric field is induced as 
expressed by the vector equation 
E = "iTx B 

where E is the induced electric potential, B is the magnetic field, and 
"u* is the water velocity. 

The induced electric potential is sensed by two electrodes, in 
contact with the water, oriented in the same direction as the induced 
field. This produces a dc signal proportional to the magnitude of the 
velocity component perpendicular to the electrode's axis. Orthogonal 
components of the velocity are measured using two pairs of electrodes 
placed with axis perpendicular to each other. 

The measureable range of velocities is from 0 to 1.5 m/sec with 
a maximum output error of one percent of full scale reading. The in- 
strument has two electrical time constant settings, 0.3 or 1.0 seconds. 

The output rms noise level is a function of the electrical time constant 
and is given by the expression 



16 



% 






1 




17 



Electromagnetic Water Current Meter EPCO 6130. 



- 1/2 

Random noise (rms) = 0.015 T (feet/sec) 

c 

where T is the electrical time constant. This relationship was validated 
c 

during calibration. 

The current meter is calibrated under steady flow conditions by the 
manufacturer but different characteristics are expected under unsteady 

flow conditions due to a change in the boundary layer on the face of the 

1 

transducer. The probe was recalibrated by oscillating it in a water tank. 

The set-up is shown in Figure 2 . The carriage on which the flow meter 

probe was mounted travels back and forth on rails . It is driven by an 

electric motor with a constant throw arm and variable rotation velocity. 

The peak carriage velocity was calculated from the tangential velocity 

of the motor arm. The ratio of carriage velocity, V^, to the velocity 

measured by the instrument, V^, was found for different angular velocities. 

This was done for each channel of the current meter by orienting each 

pair of electrodes perpendicular to the direction of flow respectively. 

The results are shown in Figures 3 and 4. 

The ratio of V /V is also plotted as a function of the frequency 
m t 

of oscillation in Figures 5 and 6. The range of frequencies measured 
ranges from 0 to 0.5 Hz; the latter being the highest frequency which 
the carriage was capable of without becoming unstable. It is noticed 
that there is a decrease in the velocity ratio with increasing frequency. 

A measured response time, or time constant, can be calculated from 
the measurements. Using the usual definition, the response time 
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Assembly for Dynamic Calibration of Water Current Meters. 

Figure 2 



MEASURED VEU 




Measured and Actual Velocities. Current Meter Serial 
Number 63 7 (Used at 69 cm Depth) 

Figure 3 
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MEASURED VEL. 




Measured and Actual Velocities. Current Meter Serial 
Number 638. (Used at 39 cm Depth) 

Figure 4 
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RATIO* MEASURED VEL. / ACTUAL VEL. 




Frequency Response of Water Current Meter 
Serial Number 637 

Figure 5 
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RATIO I MEASURED VEL. / ACTUAL VEL. 




L 

Frequency Respone of Water Current Meter 
Serial Number 638. 

Figure 6 
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a 





corresponds to the time required for the signal to decrease by 3 db from 

the true value, that is, when V /V = 0.707. Referring to Figures 5 

m t 

and 6, the instruments have an average response time of 2 seconds 
corresponding to 0 . 5 Hz . 

The probe is made of fiberglass material with dimensions 11 inches 
in length and 3/4 inches in diameter. The system senses variations in 
velocities over an area of approximately two to three probe radii from 
the transducer. This limits the spatial resolution to disturbances with 
a wave length of 11 cm, corresponding to a deep water wave period of 
0.27 seconds . 

B . WAVE GAGE 

The water surface elevation was measured with an Interstate 
Electronics Corporation SDP 201 differential pressure sensor, and a 
DP 200 wave recorder, shown in Figure 7. The pressure sensor is a 
small, unbonded strain gauge bridge. Direct current exitation for the 
bridge is supplied by a voltage regulator located in the transducer 
housing. The sea pressure is coupled by a neoprene diaphragm to a 
silicone fluid filling the interior. One part of the transducer is exposed 
directly to the interior fluid, and the other is connected to a chamber 
which is connected to the interior fluid by a length of capillary tubing. 

This arrangement acts as a hydraulic filter developing a reference 
pressure which is an average value of the external sea pressure. By 
the action of this filter, the transducer senses rapid pressure fluctuations 
only and slow changes such as tides are lost through the hydraulic filter. 
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The dc signal voltage output has a range of ^ 2 .5 volts. In the 
wave recorder, the signal is transferred to a chart paper. The recorder 
also has an outlet for an external recording system. Maximum dynamic 
range of the wave meter is + 20 feet with a linearity of one percent. 

In addition, a stilling well type mean water level indicator was 
used. This was a two inch diameter clear plastic tubing capped on 
both ends with a 1/16 inch diameter hole in the bottom and a 3/8 inch 
hole in the top for an air vent. The mean water level was measured by 
viewing the water level using a surveyor's transit from the shore and 
comparing it to a graduation on the side of the tubing. 
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lEC DP200 Portable Wave Recorder and SDP201 Wave Gage. 



IV. PRESENTATION OF DATA 



A. MEASUREMENT TECHNIQUES 

The experiments were performed at Del Monte Beach in front of 
the Naval Postgraduate School beach laboratory, on August 1, 1972. The 
waves at this location are generally highly refracted and directionally 
filtered, and break almost parallel to the shore line. At the time of 
making the measurements the breaking waves generally were of the spilling 
type with an occasional plunging type. 

It is necessary to mount the instruments on a stable platform that 
does not vibrate and can stand the forces of the breaking waves . A 
tower was built for this purpose and arranged as shown schematically 
in Figure 8. The tower was constructed of 2 inch diameter steel pipe 
with a two foot diameter flat circular base to keep it from sinking into 
the sand. It was fastened to the ground by four stays tied to screw-type 
anchors buried approximately 2 feet in the sand. The tower was tested 
for endurance before the actual taking of data for a period of 24 hours. 

The entire structure proved to be very stable; the vibration was very 
small and was neglected in the calculations. 

The transducers were arranged such that they were aligned 
vertically in order to get measurements in the same column of water. 

The two flow meters were oriented to measure vertical and horizontal 

(inshore-offshore) water particle velocities. They were mounted at 

the end of a horizontal 3/4 inch diameter pipe that extended approximately 
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Schematics of Tower Used in Taking the Measurements Showing the 
Location of The Instruments. 



equal lengths to both sides of the tower; this was for the purpose of 
minimizing any torque caused by the forces of the waves that might 
cause movement or vibrations of the probes. The wave gage was fastened 
to the base of the tower directly below the flow meter transducers. 

Figure 9 is a picture of the tower, with instruments mounted on it, 
at low tide. Measurements were made at high tide. Instruments were 
placed and recovered at low tide. Figure 10 shows the beach profile 
with the location of the tower, still water level, and the approximate 
wave breaking point at the time of recording data . 

Measurements include four velocity records (horizontal and vertical 
for each of two current meters) , sea surface elevation in the same column 
of water, and the sea surface elevation approximately 800 m offshore. A 
record of 15 minutes from the total recorded data was chosen for analysis. 
This record was taken 1 August 1972 commencing at 1705. The mean 
depth of water during the run was 1.3 meters. 

B. DATA PRE-PROCESSING 

The data was simultaneously recorded using a Sangamo Model 
3500, 14 channels FM magnetic tape recorder. The information being 
recorded was monitored continuously with a strip chart recorder and a 
volt-meter. 

The data was transcribed to an eight channel Clevit Brush strip 
chart recorder for verification using the CS-500 analog computer, and 
digitized and recorded onto a 7-track tape using anXDS-9300 computer. 
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31 



Beach Profile with Location of Measurements Site. 

V 

Figure 10 



A sample of the records from one of the flow meters and from the two 
wave recorders is shown in Figure 11. Digitization was made at a rate 
of 10 samples per second, corresponding to a sampling interval of 0.1 
second. The digitized records were transferred to a data cell in the 
IBM 360 computer for analysis. 
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Strip Chart Record Sample 
Figure 11 
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V. ANALYSIS 



A. PROBABILITY DENSITY FUNCTIONS 

The probability density functions for each of the time series was 
calculated and plotted. Comparison was made with a standard Gaussian 
distribution. The variance, standard deviation, and skewness of the 
distributions were. determined. The curves were normalized (units on 
the horizontal axis are measures of standard deviation, s) . 

1 . Sea Surface Elevation 

The pdf for the offshore wave record matches fairly well the 
normal distribution (Figure 12). This agrees with theory in that the sea 
surface elevation in relatively deep water can be approximated by a 
Gaussian distribution. 

On the other hand, the breaking wave pdf (Figure 13) was 
found to be positively skewed (skewness = 0.58). This means that 
positive values are larger but less frequent than negative values which 
are more frequent but smaller. In other words, the crests are steeper 
and more peaked, whereas the troughs are elongated and flatter. A 
schematic of such a surface is shown in Figure 14. This is how waves 
in shallow water theoretically appear. 

2 . Vertical and Horizontal Water Particle Velocities 

In relatively deep water the horizontal and vertical water 
particle velocities induced by wave motion are approximately equal for 
positive and negative directions. However, as the waves move into 
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Probability Density Function of Sea Surface Eleva- 
tion Outside Surf Zone. 

Figure 12 
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Probability Density Function of Sea Surface Eleva- 
tion Inside Surf Zone. 

Figure 13 
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shallower water they become asymmetrical and induce asymmetry in the 



velocity of the water particles , the inshore and upward velocities being 
greater in magnitude and of shorter duration than the offshore and down- 
ward velocities . This asymmetry is found in the water particle motion 
inside the surf zone and is shown by the skewness of the pdf of the 
horizontal and vertical particle velocities respectively (Figures 15 and 
16). 

B. SPECTRAL ANALYSIS 

The energy-density spectra were calculated for each variable and 
cross-spectra computed for each pair of variables inside the surf zone. 
Cross-spectral analysis was also made for the offshore and inshore wave 
records . 

The initial sampling of the records was made at a 0.1 second 
interval. However, every other sample was taken in making the compu- 
tations giving an actual sampling interval of 0.2 second and a niquist 
frequency of 2.5 hertz. Testing showed this to be a sufficiently high 
sample rate to avoid aliasing in the spectra. The maximum time lag was 
chosen to be 10 percent of the record length and the length of record 
analyzed was 15 minutes. This results in a band width resolution of 
0.0057 hertz and each spectral estimate having 20 degrees of freedom. 

1 . Offshore and Inshore Sea Surface Elevation 

The offshore and inshore sea surface elevations were 
measured using pressure wave gauges. In order to represent the pressure 



37 




cn 

XI 

Cn 

zi 

O 

u 

H 

TD 

0 

+-> 

CO 

c 

o 

w 

T3 



38 



Schematic of Sea Surface with Peaked Crests ari' 




Probability Density Function of Horizontal Water 
Particle Velocity Inside Surf Zone. 



Figure 15 
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Probability Density Function of Vertical Water 
Particle Velocity Inside Surf Zone. 

Figure 16 
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records as measures of the sea surface elevation a correction factor was 
applied to the pressure spectra. This factor is based on linear theory 
and is given by the expression 



(D(f) = 

wave 



cosh k(h - x) 
cosh kh 



(D(f) 

pressure 



where h is the bottom depth, x is the wave gauge sensor depth, k is the 
wave number, and CD (f) represents the energy -density spectra. A 
maximum limit of 100 was established for this correction factor as applied 
to the spectra. This cut-off is represented by a sudden drop in the spectra 
at higher frequencies (Figure 17). The correction factor is not applied to 
the wave spectrum in Figure 17 for the region inside the surf zone. In 
subsequent figures the correction factor has been applied to the waves 
inside the surf zone, although the validity of linear theory is highly 
questionable in this region. 

The offshore wave energy-density spectra show a maximum 
peak at a frequency of 0.19 Hz, corresponding to a wave period of 5.25 
seconds. Inside the. surf zone there is a peak at about this same 
frequency although the magnitude is smaller. Other peaks are present 
at frequencies corresponding to one-half and one-third of 0.19 Hz. 

These peaks might be caused mainly by the shoaling of the offshore 
waves. However, some kind of non-linear interaction might also be 
present as there is an increase in the magnitude of the spectra in some 
of the harmonic and sub-harmonics of the basic frequency. 
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Spectra of Sea Surface Elevation Offshore and Inside 
Breaking Zone. 



Figure 17 
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There appears to be a decrease of energy-density inside the 



surf zone compared to offshore in the frequency band of about 0.1 - 0.2 Hz. 
This decrease is apparently due to the breaking of the waves; although 
the energy-density levels for the offshore waves are questionable at 
frequencies higher than about 0.1 Hz due to the increasingly larger 
correction factor. There appears to be a shifting of energy-density to 
higher frequencies for the waves inside the surf zone. An increase 
would be expected at frequencies higher than the wave band due to 
turbulence generated during breaking. 

The two wave records show no relationship to each other in 
phase angle nor coherence. > 

2 . Sea Surface Elevation and Water Particle Velocities 

According to wave theory the wave induced horizontal 
velocity is in-phase with the sea surface elevation, and the vertical 
velocity is 90 degrees out of phase with both the sea surface elevation 
and the horizontal velocity. 

The cross -spectra of sea surface elevation with each one 
of the horizontal particle velocities at 39 and 69 cm depth were very 
similar as that shown in Figures 18 and 19. The same is true for the 
cross-spectra of sea surface elevation with the vertical velocities at 
the two depths (Figures 2 0 and 21) . 

a. Horizontal Particle Velocities 

The results show a very good coherence for each of 
the horizontal velocities in relation to the sea surface elevation in the 
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Spectra of Sea Surface Elevation and Horizontal Parti- 
cle Velocity, Depth 39 cm. 

Figure 18 
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Spectra of Sea Surface Elevation and Horizontal Parti- 
cle Velocity, Depth 69 cm. 

Figure 19 
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Spectra of Sea Surface Elevation and Vertical Parti- 
cle Velocity, Depth 39 cm. 

Figure 20 
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Spectra of Sea Surface Elevation and Vertical Parti- 
cle Velocity , Depth 69 cm. 

Figure 21 



47 



range of frequencies where most of the energy is concentrated. The co- 
herence in this range reaches a value of over 0.9 and starts decreasing 
after approximately 0.7 Hz . This is probably due to the increase in 
turbulence in higher frequencies and the general decrease in energy level. 
The phase angle shows a shift as frequency increases in the well-corre- 
lated zone of the spectrum. At low frequencies it has a value of 0 degrees, 
agreeing with linear wave theory, but increases steadily up to 90 degrees 
phase angle in the frequencies of 0.7 Hz. The phase angle becomes un- 
stable for the non-coherent region of the spectrum, 
b. Vertical Particle Velocities 

The spectral analysis of vertical velocities and sea 
surface elevation give slightly lower coherence values than for the 
horizontal velocities, but it is still over 0.75 for the range of high energy 
of the spectra. The phase angle in this case is less than the theoretically 
predicted values using linear wave theory. There is no noticeable shifting 
in phase except for frequencies higher than 0.8 Hz, where it begins to 
fluctuate, again probably due to turbulence. The area under the energy- 
density spectra is considerably less than that for the horizontal veloc- 
ities. This indicates that the kinetic energy due to vertical particle 
velocity is small. 

3 . Horizontal and Vertical Water Particle Velocities 

The cross-spectra between the two horizontal velocities 
and the two vertical velocities at the 39 cm and 69 cm depths were 
calculated with the intent of finding the relation between these parameters 
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and the distribution of kinetic energy at two different elevations in the 
water column. The two horizontal particle velocities show a very good 
coherence and 0 degrees phase angle within the significant range of the 
spectrum (Figure 22). The energy-density spectra have the same shape 
for both depths but the one calculated at depth of 69 cm has greater 
energy-density. This indicates that the kinetic energy at the shallower 
depth was smaller than at the intermediate depth. This result was not 
expected since the horizontal particle velocities according to wave 
theory are highest at the surface and decrease with depth. It is possible 
that there is a different vertical distribution of velocities after the breaking 
of the wave, but more studies are necessary in order to make any con- 
clusions . 

In the analysis of vertical velocities shown in Figure 23, a 
good coherence and 0 degrees angle of phase is also found between the 
two records in the significant energy-distribution range. The area under 
both spectra (which is proportional to the kinetic energy due to vertical 
velocities) is very small. The spectrum corresponding to the shallower 
depth is larger than the deeper one, as expected, considering that the 
vertical velocities at the bottom must be zero (if percolation is not 
considered) . 

The spectral analysis for horizontal and vertical velocity 
pairs are shown in Figures 24 and 25. The results look very similar for 
both depths. The coherence oscillates between 0.7 and 0.9 in the 
region of maximum energy. The phase angle oscillates about the value 
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Spectra of Horizontal Particle Velocity at Depth 39 cm 
and Horizontal Particle Velocity at Depth 69 cm. 

Figure 22 
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Spectra of Vertical Particle Velocity at Depth 39 cm 
and Vertical Particle Velocity at Depth 69 cm. 

Figure 23 
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Spectra of Horizontal and Vertical Velocities at 
Depth 39 cm. 



Figure 24 
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Spectra of Horizontal and Vertical Velocities at 
Depth 69 cm. 



Figure 25 
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of -90 degrees predicted by linear wave theory, but the phase shift of 
the horizontal velocity is present at both depths. The difference between 
the horizontal and vertical kinetic energy at the same depth can be 
appreciated by the difference between the areas of the energy-density 
spectra . 
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VI. CONCLUSIONS 



The values of the pdf's show that the asymmetry of the waves after 
breaking induces asymmetry in both the horizontal and vertical particle 
velocities . 

The good coherence values between the wave record and the dif- 
ferent measured velocities indicate that there is a high correlation 
between these variables; this means that the water particle velocity in- 
side the surf zone is highly wave-induced. 

The phase angle was compared to that predicted from linear wave 
theory. It shows a shifting of phase of the horizontal velocity relative 
to the sea surface elevation from 0 degree at low frequency to 90 degrees 
at higher frequencies . 

The energy-density spectra show that over 95 percent of the kinetic 
energy is due to the horizontal water particle velocity. This is in accord 
with theory for shallow water. The significant total energy is concentrated 
within a defined range of frequencies corresponding to the wave band; it 
is inferred that there is a transfer of energy to the turbulent region of the 
spectra inside the surf zone. 

From this the conclusion can be drawn that the flow motion inside 
the surf zone is not as completely disorganized as it might look, but 
rather, even after the breaking of the waves the water particle velocities 
depend on the sea surface shape and follow some determined patterns. 

More studies are necessary to get a good knowledge of the details 
particularly for different types of breaking waves . 
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